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We improve the theoretical predictions for the production of extra neutral gauge bosons at hadron 
colliders by implementing the Z' bosons in the MC@NLO generator and by computing their differ- 
ential and total cross sections in joint pr and threshold resummation. The two improved predictions 

Ch , are found to be in excellent agreement with each other for mass spectra, pt spectra, and total cross 

sections, while the PYTHIA parton and "power" shower predictions usually employed for exper- 
imental analyses show significant shortcomings both in normalization and shape. The theoretical 
1/^ ^ uncertainties from scale and parton density variations and non-perturbative effects are found to 

be 9%, 8%, and less than 5%, respectively, and thus under good control. The implementation of 
our improved predictions in terms of the new MC®NLO generator or resummed K factors in the 

(-H , analysis chains of the Tevatron and LHC experiments should be straightforward and lead to more 

f—^ I precise determinations or limits of the Z' boson masses and/or couplings. 

^. 

r-| ; I. INTRODUCTION 

fvq , Despite its impressive phenomenological success, the Standard Model (SM) of particle physics is widely believed to 

^ ' suffer from a variety of conceptual deficiencies. In particular, it provides no fundamental motivation why the strong 
On \ and electroweak interactions should be described by three different gauge groups, i.e. S'f/(3), SU{2), and U{1). Grand 
■^ Unified Theories (GUTs) allow for a unification of these groups within a simple Lie group such as SU{b), 5'O(10), 

*^ \ or Eq. Depending on the rank of the unifying group, one or several extra neutral gauge {Z') bosons appear when 
^"^ . the unification group is broken to the SM at higher scales, exhibiting the existence of additional U{\) symmetries pj. 
Anomaly cancellations and gauge invariance of quark and lepton Yukawa couplings impose a number of restrictions on 
these additional symmetries. Viable families of models, that are consistent with constraints coming from the CERN 
LEP collider, include those based on the B — L and 5*0(10) symmetries and, provided that fermion mass generation 
is not restricted to the SM Higgs mechanism and new charged fermions are allowed, also those inspired hy Eq [31 . 

If the extra Z'-bosons couple to quarks and leptons with approximately SM strength and if their mass is not too 
large, they will be produced at current and future hadron colliders and can be easily detected through their leptonic 
decay channels. The search for these particles occupies therefore an important place in the experimental programs 
C^ ' of the Fermilab Tevatron and the CERN LHC. For example, the CDF collaboration has searched the Tevatron 
Run II data for Z'-bosons in the e+e^ decay channel, using the di-electron invariant mass and angular distributions 
and setting lower mass limits of 650 to 900 GeV for a large variety of models Q. Similar constraints come from 
electroweak precision fits and di- fermion production at LEP2 [j|. Within the ATLAS collaboration, the discovery 
reach in Z' -^ e+e" decays has recently been analyzed [5] for the four classes of models defined in Ref. [^. The CMS 
collaboration has claimed a discovery reach of masses between 3.4 and 4.3 TeV for the Z' -^ /^^/i~ decay channel and 
an integrated luminosity of 100 fb~^ [f|. 

The currently available simulations for the LHC experiments rely completely on the PYTHIA Monte Carlo (MC) 
generator [7|, which is based on leading order (LO) QCD matrix elements, parton showers, and the Lund string 
hadronization model. It includes the full interference structure of new Z'-bosons with Drell-Yan photon and SM 
Z-boson exchange, and the above-mentioned phenomenologically viable models can easily be implemented [8|. The 
description of the transverse-momentum (pt) spectrum of the produced vector-boson can be improved by matching 
the parton shower to the hard emission of an extra parton 9] , but the overall normalization of the theoretical cross 
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section remains subject to large higher-order corrections and scale uncertainties. The CDF collaboration has therefore 
chosen to renormalize the generated Monte Carlo events with a correction (K) factor in each invariant mass bin to the 
next-to-next-to- leading order (NNLO) cross section [10]. However, this procedure does not lead to a correct description 
of the transverse- momentum spectrum. Note also that in principle the LO cross section no longer factorizes at NNLO, 
i.e. the i^-factors for DrcU-Yan and Z' production need no longer to be equal [2|. 

Here, we report on the implementation of extra Z'-bosons in the next-to- leading order (NLO) Monte Carlo generator 
MC@NLO fllf| , allowing to match the complete NLO matrix elements with the parton shower and cluster hadronization 
model of the Monte Carlo generator HERWIG [12]. Since the LO cross section still factorizes completely at NLO, this 
requires the implementation of the Z'-boson mass, decay width, propagator, and couplings in MC@NLO, together with 
the full interference with Drell-Yan photon and SM Z-boson exchanges. As was the case for PYTHIA, the emission of 
one additional hard parton has previously been matched to the HERWIG parton shower, albeit for SM vector-boson 
production only 13]. As an alternative, the CKKW formalism 14] for the matching of hard real emissions has been 
implemented in both PYTHIA and HERWIG for SM vector-boson production with the result that the three different 
matching formalisms were found to vary systematically in a significant way and to depend in addition strongly on the 
matching scale [iS"]. Similar results were obtained more recently in 16] for M^-boson plus jet production, where in 
addition various implementations of the MLM [l7| prescription for matching multiparton final states to parton showers 
have been compared to the methods discussed above. Here, we perform a similar study of systematic uncertainties 
by comparing the invariant-mass and transverse-momentum distributions as predicted with our implementation of 
Z'-bosons in MC@NLO with the PYTHIA parton-shower and matrix-element corrections. We also confront these MC 
predictions with a new computation of Z'-boson production in the framework of joint resummation [l8|. In addition, 
we compare the dependence of the various predictions on the unphysical renormalization and factorization scales as 
well as on the employed parton densities. The impact of hadronization corrections, dominant electroweak corrections, 
and non-perturbative effects are also studied. 

The remainder of this work is organized as follows: in SecITTlwe first describe the PYTHIA framework for Z'-boson 
production and the associated matching of parton showers with the emission of an additional hard parton. We then 
discuss our implementation of Z'-boson production in the MC@NLO generator and recall its matching procedure 
of NLO matrix elements with the HERWIG mechanism of parton showers. We also present briefly our calculation 
of Z'-boson production using joint resummation. In Sec. IIIIl we define our choice of electroweak SM parameters 
and define the parameters of our exemplary Z' model. We also discuss the various corrections that we apply to the 
production cross section, i.e. dominant electroweak corrections, next-to-leading order QCD matrix elements, parton 
showers, and hadronization. We then compare the numerical impact of these corrections and study the remaining 
theoretical uncertainties, coming from the choice of renormalization and factorization scales, the parameterization of 
parton densities, and non-perturbative effects. Our conclusions are given in Sec. IIVI 

II. THEORETICAL SETUP 

In 1984, Green and Schwarz showed that ten-dimensional string theories with Es x Eg or SO{32) gauge symmetry 
are anomaly-free and thus potentially finite [13] • Among these two gauge groups, only the former contains chiral 
fermions as they exist in the SM. After compactification, it leads to the Eg symmetry as an effective GUT group, 
that can be broken further to [20] 

Ee -> 50(10) X C/(l)^ -^ SU{5) x [/(l)^ x U{1)^. (1) 

While the Z'-bosons corresponding to the additional t/(l)^ and f7(l)x symmetries can in general mix with each other, 

Z'{0) = Z^ cos + Z^ sin 9 and Z"{e) = Z^, sin - Z^ cos 9, (2) 

we consider in this work only a TeV-scale Z;^-boson as an exemplary case, i.e. 9 = 90° in the convention of Ref. [2l| . 
and assume the Z" = Z^ to acquire its mass at considerably higher scales, as it is naturally the case in the hierarchy 
of symmetry breaking of Eq. ([1]). We will furthermore assume that 50(10) breaks down to SU{5) x U{1)^ at the 
same scale where SU{5) breaks down to the SM group SU{3)c' x SU{2)l x U{1)y with gauge couplings gs, g, and 
g'. The [/(l)x coupling 

is then directly related to the coupling g' of the weak hypercharge U{1)y by the usual group-theoretical factor -\/5/3. 
Note that as the 5O(10)-breaking scale increases from the 5[/(5)-breaking scale to the Eq or Planck scale, g^ starts 



TABLE I: Vector and axial-vector couplings of down- and up-type quarks as well as charged leptons and neutrinos in the 
-Bs-inspired Z-^ model in the convention of PYTHIA. sw is the sine of the electroweak mixing angle Qw 
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to deviate from \/5/3g' to roughly \/2/3 times this value [22] ■ Using Eq. ([3]), we can thus express g^ in terms of 
a group-theoretical factor, the SU{2)l gauge coupling g or alternatively the electromagnetic fine structure constant 
a = e^/(47r), and the squared sine of the angle dw, describing the mixing of the neutral W°- and _B-bosons to the 
massless photon and the SM Z-boson of mass mz after the breaking of SU{2)l x U{1)y to U{l)cm.- While the photon 
is protected from further mixing by the exact electromagnetic symmetry, the SM Z-boson can in general mix with 
the additional Z'-boson with an angle 4>. Its squared tangent 



tan^(j) = -^ ^ (4) 

7712 — m^ 



depends then on the eigenvalues of the Z — Z' mass matrix toi^2, which in turn depend on the vacuum expectation 
values wio and vgy^ of the S'O(IO)- and SM-breaking Higgs fields. This mixing can in general induce a couphng of the 
Z'-boson to the SM W^='=-bosons, even though they belong to different gauge groups. However, the ratio fio/^SM is 
usually large and m^ — mz ^ m2 — rnz' , so that Z — Z' mixing will be neglected in the following. The DELPHI 
collaboration has constrained (j> to be smaller than a few mrad (1.7 mrad for mz' — 440 GeV in the x-niodel) |23[. 

A. Z' production in PYTHIA 

In PYTHIA 6.403 Q, the production of extra neutral gauge bosons in hadron collisions has been implemented 
including the full interference structure with SM photon and Z-boson exchanges. The Lagrangian describing the 
interaction of the extra neutral gauge boson Z' with fermions / 

"- ^ Il^{vf-an^)IZ' (5) 
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has been expressed in terms of generalized vector [vf) and axial- vector couplings (ay), which depend in general on 
the unifying gauge group and the Higgs representations employed to break this group down to the SM gauge group. 
For the additional U{1)^ symmetry, that we use as our standard example, these couplings are given in Tab. [J for the 
down- and up- type quarks as well as for the charged leptons and neutrinos. We will be interested in the DrcU-Yan 
like production of electron-positron pairs 

qq ^ {j,Z,Z') -> e"e+, (6) 

i.e. the relevant couplings are primarily those of the five light quark flavours q = u,d,s,c,b (with masses rriq much 
smaller than the partonic centre-of-mass energy ^/s) and positrons/electrons e and to a lesser extent those of the 
other fermions, which contribute to the total decay width appearing in the Z' propagator. 

B. Matching of parton showers with LO matrix elements in PYTHIA 

The only part of the Drell-Yan like process in Eq. ([6]) that is sensitive to QCD corrections is the quark-antiquark 
(and beyond the LO the quark-giuon) initial state. It can give rise to an initial-state parton shower, that is modeled in 
PYTHIA by starting with the hard scattering partons and then successively reconstructing the preceding branchings 
in a falling sequence of spacelike virtualities Q^. They range from a maximal value Qmax = 't^zo , that is of the order 
of the hard scattering scale, to a cut-off scale Qo — ^ GeV, that is close to a typical hadron mass. 

The scale Qmax is, however, not uniquely defined. It can in particular be increased to the hadronic centre-of-mass 
energy vS, so that the parton shower (PS) populates the full phase space. However, it must then be matched to the 
QCD matrix element (ME) describing the emission of one extra hard parton 

qq^{j,Z,Z')g (7) 



at 0{aas). The emission rate for the final (and normally hardest) q ^ qg emission must therefore be corrected by a 
factor ^] 

(dg/d^)ME _ J2^=J,z.z' [^^ +u^ + 2rwf g] + interference terms 
^99 - (7.2,2') 9l^' *)- (da/dt)ps " ^^^^^^,j,[s2 + m4]+ interference terms ' ^^ 

which always lies between one-half and one. Here, s, t and u refer to the usual Mandelstam variables of the process 
in Eq. ([7|), and nij = 0. At the same order, the crossed process 

qg^i^,Z,Z')q (9) 

can occur, in which case the correction factor is 

(dg/d^)ME _ T,i=-y,z,Z'[^^ +u^ + 2m^t] + interference terms 

Since this factor always lies between one and three, the g ^ qq splitting must be preweighted by a factor three in order 
to correctly reproduce the s-channel graph qg —^ q* —^ {-y, Z, Z') q. Besides the dominant QCD radiation, we will 
also briefly investigate the effect of the relatively suppressed QED radiation (e.g. q -^ qy) as well as the corrections 
induced by the hadronization of the additional partons. 

In addition to the transverse momentum pT generated by hard emission and/or the initial-state parton shower, an 
average intrinsic transverse-momentum (fcr) can be assigned to the shower initiator in order to take into account the 
transverse motion of quarks inside the original hadron. As the shower does not evolve below Qq = 1 GeV, this same 
value is retained in PYTHIA as the default value for (A;^). However, just like Qmax, Qo = (^t) is not uniquely defined 
and furthermore closely related to the non-perturbative regime of QCD. We therefore set (fcr) — in the following. 

PYTHIA allows in principle also for the participation of multiple parton pairs in hadronic collisions. We do, 
however, not make use of this possibility, as it has little numerical effect and its description in perturbative QCD 
remains controversial. 

C. Implementation of Z'-bosons in MC@NLO 

In MC@NLO [ll|], the implementation of SM Z-boson interactions with fermions / is based on the Lagrangian [2J| 

-^/V(a/ + &/75)/^M- (11) 

cos ay/ 

For photon interactions, a/ = ey sin ^vi^ cos ^vy , where e/ is the fractional fermion charge, and 6/ = 0. The vector 
coupling a/ and the axial vector coupling 6/ are related to those defined in PYTHIA (sec Eq. (O) by a/ ^ "^f 1^ ^'^d 
hf -^ —af /A. They are combined to form the coefficients Af = ai + Ipi- and Bf — 2a/6/, that appear in the squared 
matrix elements 

\M^{qqOTqg^y,Z^e-e++X) = 1 e^ cJ Jj^^l};" 
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that have been averaged/summed over initial/final spins and colours. For the LO Drell-Yan (DY) process in Eq. ([6]), 
the colour factor is Cdy = Nc/Nl; = 1/3 and 

Tt^^\a:X = 8 t^'^? (*DY + 4y) - BiB, {tly - uly)] , (13) 

where M'^ = sdy is the invariant mass of the lepton pair and sdy , ^dy , and udy refer to the Mandelstam variables of 
the DY-process in Eq. ([6]) . The NLO QCD-corrcctions to the DY-process have been implemented in MC@NLO using 
the same convention 1241. 



One major new and technical aspect of our work is the implementation of Z'-boson interactions and exchanges 
in the framework described above. To this end, we have defined the mass and width of the Z'-boson as well as its 
couplings in the convention of PYTHIA. The squared matrix element has also been modified, 
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It includes now the squared Z'-boson exchange as well as its interferences with the photon and SM Z-boson exchanges. 
Note that it is not sufficient to change only the Z-boson mass in the existing MC@NLO implementation, but that 
it is also necessary to change the Z-boson width as well as the couplings. As a consequence, all observables change: 
the mass spectrum due to the modified width, the forward-backward asymmetries due to the modified couplings, etc. 
Furthermore, even if one changes all of these parameters, one still has 7 — Z' interference, but mo Z — Z' interference. 
We have checked numerically that these modifications induce large differences already in the total cross section. 

D. Matching of parton showers with NLO matrix elements in MC@NLO 

In HERWIG [I4I, parton showers are generated by a coherent branching algorithm with parton splittings i —^ jk, 
whose energy fractions Zj — Ej/Ei are distributed according to the LO DGLAP splitting functions. Phase space 
is restricted to an angular ordered region, which automatically takes infrared singularities into account, i.e. at each 
branching, the angle between the two emitted partons is smaller than that of the previous branching. The emission 
angles ^jk — {pj.pk)/{EjEk) — ^^L are distributed according to Sudakov form factors, that sum virtual corrections 
and unresolved real emissions and normalize the branching distributions to give the probabilistic interpretation needed 
for a MC simulation. For initial-state radiation, the parton shower follows, of course, a backward evolution. It is 
terminated when ^jk < Qg/Ef, where the space-like cut-off scale Qq is set to 2.5 GeV by default. Below this scale, a 
non-perturbative stage is imposed. In particular, a splitting of non-valence partons is enforced to allow for a smooth 
transition to the valence partons inside the outer hadron. Although HERWIG also allows for the simulation of a soft 
underlying event, we have not made use of this possibility. Since the parton shower is supposed to describe only the 
soft/collinear region, any initial-state emission outside this region, i.e. where ^ > z^, is suppressed. 

As was the case for PYTHIA, hard matrix element corrections for SM vector boson production have been imple- 
mented in HERWIG [lj| , where radiation in the region ^ > z^ can now be allowed. It is then distributed according to 
the matrix element describing the emission of an additional parton. In the soft region ^ < z^ already populated by the 
parton shower, the emission of the hardest (largest pr) parton generated so far is reweighted in order to avoid double 
counting. Note that this need not be the first emission, since angular ordering does not necessarily imply ordering in 
transverse momentum. Note also that the normalization of the total cross section |s still accurate to LO only. 

A NLO accuracy of the total cross section has been achieved in MCQNLO 11!] through an implementation of 
NLO cross sections, matched to the HERWIG parton shower. Instead of the LO matrix element implemented in 
the standard version of HERWIG, two separate samples of Born-like or standard (S) and hard emission (H) events 

are generated, that can have weight wl '' = ±1 and are statistically distributed according to the positive-definite 
standard and hard contributions to the NLO cross section Js,h- These are made separately finite by adding and 



subtracting the NLO part of the expanded Sudakov form factor and are explicitly defined in Ref. [ll| . The total cross 
section is then given by crtot = Sil'i' "^i ' (Jh + Js) /Ntot- 



E. Joint transverse-momentum and threshold resummation for Z'-bosons 



The LO matrix element predictions for Z'-production are affected by fixed order (F.O.) QCD corrections, that are 
fogarithmically enhanced, when the Z'-boson is produced close to the partonic threshold, i.e. z = M"^ j s — > 1, or 
when its transverse momentum is small, i.e. px ~* 0. These corrections must then be resummed (res) to all orders, 
which is most easily achieved in Mellin (iV) and impact parameter (&) space, and matched to the F.O. prediction by 
subtracting from their sum the perturbatively expanded (exp) resummed prediction, i.e. 
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In this way, a uniform precision is obtained, and the large theoretical (renormalization and factorization scale) uncer- 
tainties are reduced. The NLO cross section for Drcll-Yan processes (with r = M'^/S) 
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is well-known ^], and the necessary modifications for implementing Z'-bosons in the 0{al) partonic cross sections 
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have already been discussed above. 

Since the px- [2a] and threshold-enhanced contributions [2g| are both due to soft-gluon emission in the initial state, 
they may be resummed at the same time [l^, [23, UM ■ The logarithms are then organized by the function 



x(b,N) 



N 



l + r]b/N 



with b = bMe'"'/2 and N = Ne'' 
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whose form with r; = 1/4 is constrained by the requirement that the leading and next-to- leading logarithms in b and 
N are correctly reproduced in the limits 6 — > oo and N ^- oo, respectively. The choice of Eq. p^ with rj ^ 1/4 avoids 
the introduction of sizeable subleading terms into perturbative expansions of the resummed cross section at a given 
order in as, which are not present in fixed-order calculations [29|. The resummed cross section 
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can then be factorized into a regular part with 



where in the Drell-Yan resummation scheme 
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and a perturbatively calculable eikonal factor 
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on the logarithm A = Pq/t: as{iiR) In^. The anomalous dimensions T^/l^iN) are the iV-moments of the 0{as) Altarelli- 



Parisi splitting functions. Up to next-to-leading logarithmic order, the coefficients needed in ^c are 
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The usual coefficients of the QCD /9-function are 
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the number of effectively massless quark flavours is Nf, and Cp = 4/3, Ca — 3, and Tr — 1/2 are the usual QCD 
colour factors. Re-expanding the resummed cross section leads to 
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which matches precisely with the NLO cross section. 



III. NUMERICAL COMPARISON 

In this section, we compare the three different theoretical approaches to Z'-boson production at the LHC discussed 
above, i.e. the matching of LO matrix elements with parton showers as implemented in PYTHIA, the matching of NLO 
matrix elements with parton showers in our modified version of MCQNLO, and the matching of NLO matrix elements 
with our improved formalism of joint resummation at next-to-leading logarithmic (NLL) order. For a comparison of 
joint resummation with transverse-momentum and threshold resummation in Drell-Yan like processes we refer the 
reader to Ref. [H]. 

We first fix our choice of parton densities and of the strong, clectroweak, and extra gauge-boson parameters and 
demonstrate that the three theoretical predictions coincide at the LO partonic level for both the mass- and the pT- 
spectrum. Next, we show the impact of the dominant electroweak corrections by running the fixed electromagnetic 
coupling from zero momentum transfer to the mass of the Z'-boson. We then analyze separately the different levels of 
improvement made possible in the three theoretical approaches, before comparing the three "best versions" directly 
with each other. Finally, we discuss the impact of the remaining theoretical uncertainties, coming from variations of 
the renormalization and factorization scales and the parton densities, on the theoretical predictions. 



A. Choice of strong, electroweak, and extra gauge-boson parameters 



We will make predictions for pp collisions at the LHC at a centre-of-mass energy of v 5* — 14 TeV, choosing the 
CTEQ6L (LO) and CTEQ6M (NLO MS) [S^ sets as our default for the parton densities in the protons for LO 
and NLO/NLL calculations, respectively. The strong coupling constant as{fJ.R) is always computed with two- loop 
accuracy, A^-iT — 226 MeV, and identifying the renormalization scale fiR (as well as the factorization scale fip) with 
the invariant mass of the lepton pair M . As is customary, the theoretical uncertainty is estimated by varying the 
scales by a factor of two around the central value. 

For the electroweak mass, width, and coupling parameters, we use the values of the 2002 Review of the Particle Data 
Group [m, i.e. mz = 91.188 GeV, Tz = 2.4952 GeV, a = 1/137.04, and sin^ 9w = 0.23113, which are (still) used as 
default in the Z' analysis of the ATLAS collaboration. The only value that has changed in 2006 is sin^ 9w — 0.23122 
[32l |. but the numerical impact of this change remains visibly small. Using these parameters, a Z'-boson mass of 
1 TeV, running the fine-structure constant to a(l TeV)=l/124.43, and including the NLO QCD correction factor 
1 -I- as{fJ,R)/7T for Z'-decays into quarks, we compute the total width of the Z'-boson in the x-model within PYTHIA 
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FIG. 1: The LO mass-spectrum for Z' production with fixed a in PYTIfIA (triangles), MCQNLO (stars) and at parton level 
(fuU hne and circle), compared to the SM Drell-Yan background in MC@NLO (crosses). 



(see Sec. Ill A) ) to be F = 12.04 GeV. Its branching ratio to electron-positron pairs, representing the signal, is found to 
be 5.98%. These values are then passed as parameters to our modified MC@NLO and new resummation programs. 

In order to set a common theoretical basis, we show in Fig. [T] the invariant mass-spectrum of the electron-positron 
pair in LO QCD, i.e. without any parton shower, matrix element correction, resummation, or hadronization effects. 
The PYTHIA (triangles), MC@NLO (stars), and parton level (full line) mass spectra, shown in the mass range of 900 
to 1200 GeV around the mass peak of the Z'-boson at 1 TeV, coincide perfectly. For comparison, we also show the 
differential cross section for the SM Drell-Yan process in MC@NLO (crosses), which coincides with the corresponding 
LO and PYTHIA predictions and represents the dominant (irreducible) background to the Z'-boson signal. Far from 
the resonance region, the two mass spectra would, of course, coincide. The transverse momentum of the lepton pair 
(pt) is exactly zero in PYTHIA and at parton level, while the forced splitting of non-valence partons in MCQNLO 
induces a distribution that extends to non-zero values of px even when the parton shower is switched off. Since 
this distribution is unphysical, we do not show it here. When integrated over px, the MC@NLO total cross section 
coincides, however, with PYTHIA's and the one at parton level. 

The dominant electroweak corrections can be resummed by running the fixed value of the fine-structure constant 
a = 1/137.04 in Fig. [T] to the value a — 1/124.43 at the mass of the Z'-boson in the MS-scheme. The ratio 
a^(M)/a^(0) then induces an increase of about 22% in the cross sections shown in Fig. [TJ In the following, we will 
always use a running value of a. 



B. Numerical results with PYTHIA 



In Fig. [21 we show the three different ways of improving on the parton-level predictions that are implemented in 
PYTHIA, where one can either add QCD parton showers in the soft and collinear regions (circles), corresponding 
to the leading- logarithmic (LL) approximation, or in the full phase space (triangles), which however overestimates 
the cross section, so that it must be renormalized to the matrix element describing the emission of an additional 
hard parton (squares). By definition, neither the total cross section nor the mass spectrum are changed and remain 
accurate to LO only, as can be seen from the fact that all histograms on the left-hand side of Fig. [5] coincide, after 
normalization to the LO prediction and within the statistical error bars, with unity. The right-hand side of Fig. [2] 
demonstrates that the zero transverse momentum of the .Z'-boson is smeared by the parton showers. It then peaks 
around 3 GeV and extends up to 700 GeV (see insert) and even beyond, if the available phase space is opened to the 
full hadronic centre-of-mass energy as proposed by the "power shower" prescription. Note that strictly speaking the 
LL approximation is, of course, no longer valid there. Normalizing the "power shower" to the correct QCD matrix 
element describing one real emission brings the integrated cross section back into agreement with the LO prediction. 
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FIG. 2: Mass (left) and transverse-momentum spectra (right) with PYTHIA with soft/colhnear QCD parton showers (circles), 
QCD parton showers populating the full phase space (triangles), and after adding LO matrix element corrections (squares). 
The mass spectra have been normalized to the LO QCD prediction. 

as the increase of the cross section at large pT is compensated by a reduction at small pr- These predictions may be 
improved further in PYTHIA by adding QED parton showers, which are, however, formally suppressed by a factor 
of a /as, so that their influence would not be visible in Fig. [2] The emitted additional partons may furthermore 
hadronize, which is modeled in PYTHIA with the Lund string model, but we have checked again that the shapes of 
the distributions in Fig. [2] would not change significantly. 

C. Numerical results with MCONLO 

With MC@NLO it is possible to correct not only the transverse-momentum shape, but also the normalization 
of the total cross section. In our implementation, this is now also possible for Z' production at the LHC. This is 
demonstrated in Fig. [3l where we show on the left-hand side that the NLO correction factor 



K 



dcTNLO 
dcTLO 



1.26 



(29) 



at the resonance (M — 1 TeV) is quite significant and depends also slightly on the invariant mass of the lepton pair. 
By definition, the normalization remains again unchanged by the HERWIG parton shower, which affects, however, 
strongly the pT^-spectrum (right). While the fixed-order prediction diverges as px — > 0, its logarithmic singularity is 
effectively resummed by the parton shower and leads to a smooth turnover with a maximum at around 8 GeV, i.e. at 
a value that is considerably larger than in the case of PYTHIA. Note that the parton shower in MC@NLO replaces 
part of the NLO contribution, so that switching it off (not shown) does not lead to a fully correct NLO prediction. 
As was already the case in PYTHIA, we have checked that adding QED parton showers or hadronization, which is 
modeled in HERWIG with the cluster model, docs not lead to visible changes in the distributions of Fig. [31 



D. Numerical results with joint resummation 



Theoretically, the most precise predictions are obtained if logarithmically enhanced soft and coUinear parton emis- 
sion is analytically resummed, both close to threshold and close to pr — 0, as described in Sec. Ill El above. As 
can be seen from Fig. [5] (left), the NLO iiT-factor (dashed) is increased further by the resummed contributions, in 
particular in pure threshold resummation (dotted), which is only approximately described by joint resummation, as 
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FIG. 3: Mass (left) and transverse-momentum spectra (right) after matching the NLO QCD corrections to the HERWIG QCD 
parton shower (triangles) . The mass spectra have been normalized to the LO QCD prediction. 
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and the mass spectra have been normalized to the LO QCD prediction. 



we are still relatively far from the production threshold at ^/S — 14 TeV. Only for larger values of M would the joint 
and threshold resummed predictions coincide. As the additional increase in K is only of the order of a few percent, 
the prediction is nicely stabilized. We have checked that the iiT-factors of 1.26 at NLO and 1.28 at NNLO [iO| for 
Drell-Yan lepton pairs with mass 1 TeV coincide precisely with our NLO and threshold resummed results. Fig. [4] 
(right) demonstrates that the true NLO py-distribution (dashed) diverges as pt —* 0, as it must, and does indeed not 
exhibit the unphysical maximum of MC@NLO without parton showers around 10 GeV. Resummation leads again to 
a smooth turnover with a maximum at around 8 GeV. Note that the jointly resummed prediction (full) follows the 
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FIG. 5: Mass (left) and transverse-momentum spectra (right) in PYTHIA with LO matrix elements matched to QCD parton 
showers (circles), in MC@NLO with NLO matrix elements matched to the HERWIG QCD parton shower (stars), and after 
matching the NLO QCD corrections to joint resummation (full line) . The mass spectra have been normalized to the LO QCD 
prediction, and the renormalization and factorization scale uncertainties in the resummed predictions are indicated as shaded 
bands. 

one with pr-resummation (dotted) very closely, as is to be expected from the definition of the variable x ^nd the form 
of the function Qc (see above) . 



E. Comparison of numerical results and theoretical uncertainties 



We now confront the "best versions" of the three different theoretical approaches with each other, superimposing 
in Fig. [S] the LO matrix element corrected predictions obtained with PYTHIA's "power shower" (circles), the NLO 
matrix element corrected predictions obtained with MC@NLO's parton shower (stars), and the jointly resummed 
prediction matched to the NLO matrix elements (full line). The correction factors for the mass spectra (left), which 
have been normalized to the LO QCD prediction, show only a very weak mass dependence. We have multiplied the 
PYTHIA mass spectrum by hand with a global iiT-factor of 1.26. Otherwise, within the statistical error bars the 
PYTHIA K-iactoT would just be unity, since the normalization of the total cross section is changed neither by the 
parton shower nor by the LO matrix element correction, that serves to bring the "power shower" back into agreement 
with the LO QCD prediction. The MC@NLO i^-factor agrees almost perfectly with the one of joint resummation, 
since we saw in Fig. S] (left) that threshold and joint resummation lead only to a very modest increase of the NLO 
ii'- factor. The theoretical uncertainty induced in the resummed prediction through the simultaneous variation of the 
renormalization and factorization scale by a factor of two around the central scale M is also considerably smaller 
than the K-iactoi, indicating a nice stabilization of the theoretical prediction. The p^-spectra (right) are, for all 
three improvements, no longer divergent. The PYTHIA prediction rises and falls rather steeply around its maximum 
at 3 GeV, whereas the MC@NLO and resummed predictions rise and fall more slowly around the peak at 8 GeV, 
which has furthermore a slightly smaller cross section. The agreement between MC@NLO and joint resummation is 
impressive, in particular for a scale choice of M/2 (upper end of the shaded band) at low pr and M at intermediate 
Pt- We can therefore conclude that our implementation of Z' production in MC@NLO reaches almost the same level 
of precision as our joint resummation calculation, but offers the additional advantage of an easy implementation in 
the analysis chains of the LHC experiments. 

The scale uncertainty of the total cross section (integrated over all transverse momenta and over the invariant mass 
in the range from 900 to 1200 GeV) is shown in Fig. E The LO QCD prediction (full) agrees with the PYTHIA 
prediction (dot-dashed) at the same order, as the total cross section and its scale dependence is not modified by the 
parton shower. Since as{fJ-R) does not enter the calculation at this order, the full scale dependence is in fact due to 



12 







PP -^ 


y, Z, Z' -^ r 


r at the LHC 






-■■ 










420 


LO 












NLO 






_ 






Joint Resummation 




400 


- 






PYTHIA 

MC at NLO 




380 








360 


_ 


" ' ' — 


— ~l:v_ 






tT 


- 








" —^ ■--. 


320 


r 










300 


-^ 


^__^ 








280 


- 




^^^^^— -^ 






260 


_ 








^^^—-^ 




" 




II 



F R 
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LO matrix elements to the PYTHIA parton shower (dot-dashed), and NLO matrix elements to the HERWIG parton shower 
(long dot-dashed). 

the factorization scale. The LO scale dependence does, however, not give a reliable estimate of the theoretical error, 
since the NLO cross section (dashed) is considerably larger. At NLO, the factorization scale dependence is reduced 
as expected, but as{fiR) makes its appearance, so that an additional renormalization scale dependence is introduced. 
The total NLO scale dependence is reduced to 9% (vertical lines), once the leading and next-to- leading logarithms 
(NLL) are resummed (dotted). The MC@NLO prediction (long dot-dashed) agrees with the one at NLO (dashed) 
for the central scale, but has a weaker scale dependence due to the resummation of leading logarithms in the parton 
shower |33| . Since both the type {pt vs. joint) and order (LL vs. NLL) of the resummed logarithms differ between 
MC@NLO and joint resummation, the two scale variations need not (and do not) coincide. 

We estimate in Fig.[7]the theoretical uncertainty coming from different parameterizations of parton densities. Since 
the invariant mass M of the lepton pair is correlated with the momentum fractions Xa.b of the partons in the external 
protons, the normalized mass spectra (left) are indicative of the different shapes of the quark and giuon densities in 
the CTEQ6M (full) and MRST 2004 NLO [H] (dashed) parameterizations. The latter also influence the transverse- 
momentum spectra, which are slightly harder for MRST 2004 NLO than for CTEQ6M. The shaded bands show the 
uncertainty induced by variations along the 20 independent directions that span the 90% confidence level of the data 
sets entering the CTEQ6 global fit [3g. It remains modest, i.e. about 8%, and is thus slightly smaller than the scale 
uncertainty of 9%. 

The uncertainty at low transverse momenta coming from non-perturbative effects in the PDFs is usually parame- 
terized with a Gaussian form factor describing the intrinsic transverse momentum of partons in the proton. We show 
in Fig. [5] the effects coming from three different choices of the form factor, i.e. those of Ladinsky-Yuan (LY-G) [33], 
Brock-Landry-Nadolsky-Yuan (BLNY) [11], and Konychev-Nadolsky (KN) (39tl, on the quantity 



A 



da("---+NP)(;i, 



flF 



M)-da'-''''--'>{fiR^fiF = M) 



dcr(™*'-)(^fl, ^ fiF ^ M) 



(30) 



It is obvious that these non-perturbative contributions are under good control, as their effect is always less than 5% 
for Pt > 5 GeV and thus considerably smaller than the scale and PDF uncertainties. 



IV. CONCLUSIONS 



In summary, we have improved the theoretical predictions for the production of extra neutral gauge bosons at 
hadron colliders by implementing the Z' bosons in the MCQNLO generator and by computing their differential 
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FIG. 8: Variation (in percent) of the transverse-momentum spectrum after matching the NLO QCD corrections to joint 
resummation with CTEQ6M parton densities for three three different choices of a non-perturbative form factor. 



and total cross sections in joint px and threshold resummation. The two improved predictions were found to be in 
excellent agreement with each other for mass spectra, pr spectra, and total cross sections, while the PYTHIA parton 
and "power" shower predictions usually employed for experimental analyses show significant shortcomings both in 
normalization and shape. The theoretical uncertainties from scale and parton density variations and non-perturbative 
effects were found to be 9%, 8%, and less than 5%, respectively, and thus under good control. The implementation 
of our improved predictions in terms of the new MC@NLO generator or resummed K factors in the analysis chains 
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of the Tevatron and LHC experiments should be straightforward and lead to more precise determinations or limits 
of the Z' boson masses and/or couplings. While we have shown numerical results for Z' bosons associated with the 
C/(l)^ gauge symmetry, our calculations are completely general and easily applicable to different grand unification, 
extra-dimensional, or other models. Our modified MC@NLO program has been endorsed by the original authors and 
can be obtained, like the joint resummation program, at http : //Ipse . in2p3 . f r/klasen/sof tware . 
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